FERMILAB-PUB-00 /346-T 
hep-ph/010102^ 



The Semileptonic Decays B — » ixlv and D —> ttIu from Lattice QCD 



Aida X. El-Khadra, 1 Andreas S. Kronfeld, 2 Paul B. Mackenzie, 2 Sinead M. Ryan 3 and 

James N. Simone 2 
1 Physics Department, University of Illinois, Urbana, IL 61801, USA 
2 Fermi National Accelerator Laboratory, Batavia, IL 60510, USA 
^School of Mathematics, Trinity College, Dublin 2, Ireland 
(January 2, 2001) 



Abstract 

We present a lattice QCD calculation of the form factors and differential decay 
rates for semileptonic decays of the heavy-light mesons B and D to the final 
state ttIu. The results are obtained with three methodological improvements 
over previous lattice calculations: a matching procedure that reduces heavy- 
quark lattice artifacts, the first study of lattice-spacing dependence, and the 
introduction of kinematic cuts to reduce model dependence. We show that the 
main systematics are controllable (within the quenched approximation) and 
outline how the calculations could be improved to aid current experiments in 
the determination of \V u b\ and IV^I. 



PACS numbers: 12.38.Gc, 13.20.He, 12.15.Hh 
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I. INTRODUCTION 



Processes involving weak decays of B and D mesons are of great interest, because they 
yield information on the more poorly known elements of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix. Semileptonic decays have traditionally been used to determine the CKM 
matrix, for example, V u d (through nuclear /3-decay), V us (Ki 3 ), V c b (B — > D^*Hu), and 
V u b (b — > ulv) 0. In the first three cases flavor symmetries (isospin, SU(3) flavor, and 
heavy quark symmetry, respectively) greatly simplify one's theoretical understanding of the 
hadronic transition matrix elements. In the symmetry limit, and at zero recoil, current con- 
servation ensures that the matrix elements are exactly normalized. Even when estimates of 
the deviations from the symmetry limit are difficult to calculate reliably, the deviations tend 
to be small. Thus, the overall theoretical uncertainty on the decay process is under con- 
trol. Given good experimental measurements, this procedure then determines the associated 
element of the CKM matrix. 

For semileptonic decays of charmed or 6-flavored mesons into light mesons there are no 
flavor symmetries to constrain the hadronic matrix elements. As a result, the errors on \V u b\ 
are currently dominated by theoretical uncertainties and are not well known |]T|. For the 
same reason the best value for \V c d\, at this time, comes from neutrino production of charm 
off of valence d quarks (with the cross section from perturbative QCD), rather than from 
the semileptonic D decays. In this paper we take a step towards reducing the theoretical 
uncertainty by using lattice QCD to calculate the form factors for the decays B — > tiIu and 
D — > nlu. Although our results are in the quenched approximation, we introduce several 
methodological improvements that carry over to full QCD. Moreover, this work is the first 
to study the lattice-spacing dependence of the form factors. 

There is a considerable ongoing experimental effort on this subject, which will lead to 
measurements of the differential decay rates. For B — > nlu, 

dT G 2 F \V ub \ 2 2m B p*\f + (E)\ 2 

dp 2471- 3 E [ ' j 

where E = p n ■ p B /m B is the energy of the pion in the rest frame of the B meson, and 



p = J E 2 — m\ is the magnitude of the corresponding three-momentum, (p^ and p B are 
four-momenta. For D — > nlu, replace V u b with V c d, m B with mp, and p B with The 
non-perturbative form factor f+{E) parametrizes the hadronic matrix element of the heavy- 
to-light transition, 



(7T( P ^\B(p B )) = f + (E) 



m\-m\ 

Pb+Ptt 5 — q 

q 2 



\fo(E)^^q», (1.2) 
q 2 



where V M is the charged b —>■ u vector current, and q = ps—p-n is the momentum transferred 
to the leptons. For reasons that are made clear below, we prefer to consider the form factors 
f + and fo as functions of E. This kinematic variable is related to the more common choice 
q 2 = m\ + m 2 — 2m B E. The contribution of fo to the decay rate is suppressed by a factor 
(mi/m B ) 2 so we shall present the rate given in Eq. ( |1 . 1| ) . In the decay B — > tttv both form 
factors are important, however, so both are tabulated below, in Sec. |VI| 
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The first determinations of \V u b\ came from the rate of the inclusive semileptonic decay 
B — > X u lv. In general, inclusive rates can be described model-independently through an 
operator product expansion (OPE), leading to a double series in AqcdA^ an d a s {m b ) [0. 
Thus, they are subject to non-perturbative and perturbative uncertainties. In particular, 
one requires the quantities A, Ai, and Ag, which are defined in the heavy-quark effective 
theory^] The huge charm background in B — > X u lv must be eliminated by imposing a cut 
either on the charged lepton energy M , on the hadronic invariant mass , or on q 2 M . Such 
cuts narrow the kinematic acceptance and may, therefore, increase sensitivity to violations 
of quark-hadron duality, which is hard to quantify 

The differential rates of exclusive decays offer an alternative route to and |V^|. 
This method is limited, however, by uncertainties in the form factors, such as f+(E) in 
Eq. ( |1.1| ). In the case of D decays, the E dependence of the rate has been measured only 
for D — > Klv @. The FOCUS collaboration || will improve that measurement and also 
should be able to measure the E dependence in the Cabibbo-suppressed mode D — > ttIu. 
First measurements of the branching ratios for B — > nlv and B — > plv have been presented 
by the CLEO collaboration 0. The form factors for all these processes are calculable with 
lattice QCD. Here we concentrate on calculating the form factors for B — > irlv (and similar D 
decays). The branching ratio is not as large as for B —>■ plv, and there are other experimental 



difficulties [10]. On the other hand, with vector mesons several form factors enter into the 
decay rate. Furthermore, one might expect greater uncertainties for the p (and to and </>) 
from the quenched approximation, because of their non-zero hadronic widths. 

With lattice QCD a very pressing issue is to understand the systematic uncertainties. 
Indeed, an important justification for using the quenched approximation is that the savings 
in computer time allow us to study the other systematic uncertainties in detail. To control 
systematic errors we apply three main methodological improvements in this paper: we nor- 
malize the heavy-quark action and current in a way that reduces heavy-quark discretization 
effects, we have three different lattice spacings to study any remaining discretization effects, 
and we introduce kinematic cuts to avoid model dependence. 

First, let us consider the discretization for the heavy quark. At the lattice spacings, a, 
currently in use, the large mass of the b quark means that m^a > 1. To control lattice 
spacing effects, we adopt the approach of Ref. |L1|, which takes an improved action for 



Wilson fermions, but adjusts the couplings in the action and the normalization of the current 
so that the leading and next-to-leading terms in the heavy-quark effective theory (HQET) 
are correct. By applying HQET directly to lattice observables, one can show that the 



heavy-light meson has small discretization effects [12], in our case of order cXs-^qcd/^q, 
Q^Aqcdo, (Aqcd/ttiq) 2 , and (Aq CD a) 2 . These normalization conditions allow us to perform 
our calculations directly at the physical mass rriQ = nib- This approach has already been 
successfully applied in calculations of B and D meson decay constants by four groups Hl3]-[16| 
and in calculations of the form factors for B — > D^*'lv at zero recoil [|I7]]. Work on B — > 
irlv by two other groups |flS| -|20fl with the same action (but different lattice currents) has 
used normalization conditions designed for light quarks, which suffer from errors of order 



1 A new method for calculating A, Ai, and A2 can be found in Ref. El. 
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a s rriQa |L8| or (mga) 2 |19|,p0 . 



To reduce these effects their calculations have been carried 

or 1.7 GeV < mp < 



19 



out with pseudoscalar meson masses 1.2 GeV < mp < 2.0 GeV 
2.6 GeV p0[| . We have not been persuaded that HQET can be used to guide the extrapolation 
from there back up to m B = 5.3 GeV. 

Second, there are cutoff effects of order a s ka and (ka) 2 from the light quark, where k 
is the momentum of the light quarks inside the mesons. For light or heavy-light hadrons 
at rest, the momentum k « Aqcd, so these effects are of the same kind as some of those 
considered above. In the semileptonic decay, however, one has a light daughter hadron with 
non-zero recoil momentum, which gives rise to lattice spacing errors with k = \p n \. To study 
this systematic error, we carry out the calculation at three different lattice spacings, and 
check the dependence of our results on a. We can then restrict our final results to small 
enough recoil momenta, so that discretization effects remain under control. Our test of the 
lattice spacing dependence is the first in a lattice calculation of semileptonic form factors. 

Third, we do not use models to extend our kinematic reach to high pion energy (i.e., 
low q 2 ), in contrast to previous work [|18| -|2"0"|| . The extrapolation would rely on the worst of 
our data: not only do discretization errors increase with p w a, but statistical errors do too. 
Therefore, we quote the differential decay rate over the range where systematic uncertainties 
from the lattice are under control. In particular, we define 



(Pmin? Pmax) 



d PP y + (E)\ 2 /E. 



;i.s) 



The upper limit is chosen to rein in the discretization and statistical uncertainties. The lower 
limit cuts out a region where extrapolations in p and light quark mass are difficult. Then, 
assuming a massless charged lepton, one can combine T B with experimental measurements 
to determine the CKM matrix via 



\V ub \ 2 



12tt 3 



1 



G 2 F m B T B (p 



Pmax dT 

dp 

Pmin U-P 



B^TT 



and, similarly, 



\v cd \ 2 



12tt 3 



G 2 F m D T D (p 

miri) Pmax J 



Pmax dT 

dp 

Pmin dp 



D — >7T 



;i.5) 



Our final result, showing the integrand of Eq. ( |1.3| ) for B and D — > tiIu, is in Fig. [I]. 
The shaded regions indicate the range of pion momentum over which we can control the 
uncertainties. Integrating over this region, we find 



T B (0A GeV, 1.0 GeV) = 0.55; 
T D (0A GeV, 0.925 GeV) = 0.23; 



-0.15 
-0.05 
-0.04 
-0.02 



-0.09 +0.09 
-0.12 -0.02 
-0.01 +0.06 
-0.05 -0.02 



± 0.06 ± 0.09 GeV 4 , 
± 0.03 ± 0.03 GeV 4 . 



;i.6) 
;i.7) 



where the first uncertainty is statistical, and following four are systematic and come from 
chiral extrapolation, lattice spacing dependence, matching to continuum QCD, and the sum 
in quadrature of several other uncertainties. The last includes an estimate of the uncertainty 
from converting lattice units to physical units, which partly reflects uncertainty from the 
quenched approximation. In addition to these uncertainties, which are quantifiable within 
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FIG. 1. The differential decay rate (without momentum- independent factors) as a function 
of p = \Ptt\, for (a) B — ► ttIv and (b) D — ► 7r^. The solid error bars show the statistical uncertainty 
and the dotted ones show the sum in quadrature of statistical and systematic uncertainties. 
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the quenched approximation, there may be an additional error from quenching as large as 
10-20 percent on T B and T D . 

At low momenta the experimental rates go to zero, so no information is lost by making the 
cut at p m i n = 0.4 GeV. For semileptonic D decays the high-momentum cut is already at the 
kinematic endpoint (m 2 D — m 2 T )/2mD = 0.925 GeV. A high-momentum cut at p max = 1-0 GeV 
is, however, an obstacle to determining |V4&|, since semileptonic B decays usually produce 
harder pions. Although the cut does reduce the overlap between our lattice calculation and 
experimental results, the results presented here are model independent (apart from quench- 
ing). As experimental and lattice results improve over the next several years, the range of 
pion momentum should widen and can be selected to optimize the combined experimental 
and theoretical uncertainty. 

This paper is organized as follows. Section [IT] contains a discussion of the lattice action 
and vector current for heavy quarks. The lattice calculation of the matrix elements is 
described in Sec. [TXJ]. Section [TV] describes an interpolation in pion three-momentum and an 
extrapolation in light quark mass, which are needed to obtain the form factors. The former 
is a special feature of these decays; it interacts with the chiral limit, and together these lead 
to the cuts given in Eqs. ( |1.6| ) and ( |1.7| ). We discuss quantitatively the systematic errors 
on Tb and T B in Sec. [V]. The analysis of B and D decays is essentially the same. Results 
for the form factors are tabulated in Sec. [VI|. Section |V11| compares our methods and results 
to previous (and ongoing) work |19|-^I|. Section |VII1| concludes. 



Preliminary results of this analysis have been presented in Refs. [p2p3| . Phenomeno- 
logical implications of D decays, especially for comparing D — > tiIu and D — > Klv as in 
Ref. |[24||, will appear in another publication. 



II. CONTINUUM AND LATTICE MATRIX ELEMENTS 



The continuum matrix element of the flavor- changing vector current, V M = uij^b, is 
parametrized by two independent form factors, for example those in Eq. ( |1.2|) . In considering 
the chiral and heavy-quark limits, it is more convenient to write the matrix element as 

(7r( P7r )\V»\B(p B )) = V2^[v»f\\(E)+if ± f ± (E)] , (2.1) 

where v = Pb/^b is the four-velocity of the B, and p± = p n — Ev is the pion momentum 
orthogonal to v. The traditional form factors /+ and fo are related to f\\ and f± by 

(2ms)- 1 / 2 



U(E) 
HE) 



\ / 2m B 



m 



B 



mi 



f\\(E) + (m B -E)f\(E) 
(m B -E)f ll (E) + (E 2 -ml)f ± (E) 



(2.2) 
(2.3) 



At q 2 = it follows from these formulae that f+ = fo, which is necessary from Eq. ( |1 . 2| ) . 

There are several good reasons to focus the numerical analysis on f» and f±. 
consideration of chiral and heavy-quark symmetry yields the expectation for m n , E 

fB^/m^ 



■/ll 



fi 



V2f« ' 

fB*y/m~B~* 

V2U 



■9bb* 



2m i 



mi 



First, 
-+ 

(2.4) 
(2.5) 
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through order 1/mq in the heavy-quark expansion [p5|| . Here fs*, and /,r are decay 



constants, and gBB*-w is the B-B*-n coupling. Although we do not use these results to 
constrain the needed chiral extrapolation of our data, they do show us that f\\ and f± 
behave differently as m n is reduced to its physical value. (Recall q 2 = m 2 B + m 2 — 2m B £.) 
Furthermore, f\\ and f± have a simple description in the heavy-quark effective theory |25fl . 
so they are natural quantities to study in the lattice method of Refs. |~T],[T2|, or when using 



lattice NRQCD [FZIJ. Finally, they emerge directly from the lattice calculation, so it is 



simpler to analyze them separately, forming the linear combinations in Eqs. (|2.2|) and ( |2.3| ) 
at the end. 

For the light quarks we use the Sheikholeslami-Wohlert (SW) action [pSfl , with the cus- 
tomary normalization conditions for m q a — > 0. The SW action has an extra coupling csw, 
sometimes called the "clover" coupling, which can be adjusted to reduce the leading lattice- 
spacing effect of Wilson fermions. In practice, we adjust csw according to tadpole- improved, 
tree-level perturbation theory [P7|1 , so the leading light-quark cutoff effect is of order a s ka. 

We also use the SW action for the heavy quark, but its two free parameters, the bare 
mass mo and clover coupling csw, are adjusted to maintain good behavior in the heavy- 
quark limit |fn}| . This goes as follows: on-shell lattice matrix elements can be described by 



a version of HQET [T2j , with effective Lagrangian (in the rest frame) 

- hD 2 h hiH-Bh , . 

£ H qet = rnihh + — 1 1 , (2.6) 

ZUI2 ZTTlQ 

where h is the heavy-quark field of HQET, and B is the chromomagnetic field. The "masses" 
mi, m 2 , and m B are short- distance coefficients; they depend on m and csw ( an d the gauge 
coupling). Fortunately, matrix elements are completely independent of vii\ ||12|| , so we ad- 
just mo and csw to tune m 2 and mg to the b (or c) quark. In practice, we tune m 2 
non-perturbatively, using the quarkonium spectra, and me with the estimate of tadpole- 
improved, tree-level perturbation theory |^7j . 

The lattice current is constructed according to the same principles. We distinguish the 
lattice current V from its continuum counterpart V and take 



V = ^Z V uuZ v »*J'fV b (2.7) 

where the rotated field [|TT|] 

= [1 + arf l7 • Aat] ip q , (2.8) 

and ip q is the lattice quark field (q = u, b) in the SW action. Here Z?i at is the symmetric, 
nearest-neighbor, covariant difference operator. In Eq. ( |2.7j ) the factors Z vqq , q = u,b, 
normalize the flavor-conserving currents. In practice, they are computed non-perturbatively. 
Matching the current V^ 1 to HQET requires further short- distance coefficients: 

V" = (# + C^qh + V ^ ± h - + . . . , (2.9) 

zm 3 

where the symbol = implies equality of matrix elements, and q is a relativistic (continuum) 
anti-quark field. At the tree level rfy 1 = 1, (y 1 = 0. Also, further dimension-four operators, 
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whose coefficients vanish at the tree level, are omitted from the right-hand side of Eq. Q2.9Q . 
This description is in complete analogy with that for the continuum current, namely, 



V M = {r] V + (v)v^qh + r] V qi^h 



2mr 



+ 



(2.10) 



Indeed, the HQET operators are the same. On the other hand, the radiative corrections 
to the short- distance coefficients in Eqs. ( [2.9|) and Eqs. ( |2.10| ) differ, because the lattice 
modifies the physics at short distances. 

one can deduce how to compensate 



By studying the form factors in HQET, as in Ref. |25 



for the mismatch between short- distance coefficients rjy 1 ^ and £y M} for the lattice and r/y 
and f° r the continuum. HQET matrix elements have form factors 



(n\qh\B)= n (E), 
(rc\qtY ± h\B)=p^±(E), 

so, leaving aside the dimension-four operator qi^p^h for now, 

f ± (E) = (r 1v + ( v )v ± (E). 
By the same reasoning, form factors calculated with the lattice current V satisfy 

/ii at (£) 



(^ at + Cy>±(£). 



(2.11) 
(2.12) 



(2.13) 
(2.14) 



(2.15) 
(2.16) 



Up to lattice artifacts of the light degrees of freedom the HQET form factors <p\\ and ip± are 
the same in Eqs. (fDl) and (gTp and in Eqs. (gTp and (grp . Thus, 



f\\(E) 
f±(E) 



Pv ]{ f{\(E), 
PvJ^(E), 



(2.17) 
(2.18) 



where pv« = w/^l/S Pv ± = (vv + Cv)/(^l/ t + Cy*)- Because these factors arise from short 
distances, in practice we compute them in perturbation theory to one loop. We find these 
short- distance corrections to be very small. 

Finally, the free parameter d\ in Eq. ( |2.8| ) can be adjusted to tune l/m 3 to BxjvriQ. 
In the present calculations, we adjust d\ with the estimate of tadpole-improved, tree-level 



perturbation theory, as explained in Ref. [ I I 



With these normalization conditions the leading term in the heavy-quark expansion is 
correctly obtained, up to neglected higher-order corrections to pv„ and pv ± - The associated 
error should be much smaller than our other uncertainties, because most of the short- distance 
normalization is handled non-perturbatively, through the factor \J Zy^ Z V bb . Similarly, the 
1/rriQ term in the heavy-quark expansion is correctly obtained, up to neglected loop cor- 
rections to csw and d\, and to dimension-four operators neglected in Eq. (|2.9|). Here the 
associated error depends on rriQa. When rriQa > 1 it is formally of order a s AQCD/^Q, but 



8 



when rriQa < 1 it is formally of order o! s AqcdO- In the work reported here, such corrections 
are smaller than, or comparable to, other uncertainties. 

In lattice QCD the required matrix elements and thence the form factors are calculated 
from correlation functions. In particular, the three-point correlation function for the B — > tt 
transition is 

Cf (p, k,t f , t„, U) = J2 e- tpx e^ k - p >y(0\O B (x, t f )%j^ u (y, t,)0t (0, t t ) |0>, (2.19) 

x,y 

where Ob and O w are interpolating operators for the B and 7r mesons. In the limit of large 
time separations, the correlation function becomes 

(p, k, tft ts , U ) = 4/^1/2 (^y)) e -^ ts - ti}e - EB{tf - ta) + . . . } (2 20) 

where Eb {E n ) is the energy of a B (tt) meson with momentum k (p). The energies and 
the external line factors Z n and Zb can be calculated from two-point correlation functions 

C®(p,t) = J2e~ ipx (0\O H (x,t)OU0,0)\0), (2.21) 

X 

where H is 7r or B, and for large \t\ one has 

C®(p,t) = Z H e- EH ^ + ■ ■ ■ . (2.22) 

By time reversal (i?(fc)| 1 l r fe 7 Ai \I/ u |7r(p)) = (n(p)\^/ u j^ b \B(k)), so in the rest of this paper 
we do not distinguish the two matrix elements. 

To summarize this section, let us review the steps needed to obtain the physical form 
factors /+ and fo. First we obtain E n and 

for several values of p, directly from fitting the lattice correlation functions to the time 
dependence given in Eqs. ( [2.20|) and Q2.22|) . The normalization factors Zy^ and Z ybb are 
computed from zero-momentum, flavor- conserving correlation functions. The radiative cor- 
rection factors pv appearing in Eqs. ( |2.17| ) and Q2.18Q are computed with perturbation theory. 
These ingredients are combined to form 



/||(E) = p v ^Z V u U Z vbb pE^F\\(p), (2.25) 



f ± (E) = p v± yJZ V u*Z vbby /2E w aF ± (p), (2.26) 

with E = E n . From the calculated values of p we then interpolate to a fiducial set of 
momenta. The form factors /y and f± are extrapolated to the physical light quark mass. 
With the light quark corresponding to strange we check also for lattice spacing effects. 
Finally, the combinations / + and /o are formed from the extrapolated f± and /11 with 
Eqs. (|2.2j) and (|2.3|) and physical meson masses. 



9 



TABLE I. Input parameters to the numerical lattice calculations, together with some elemen- 
tary output parameters. Error bars on the outputs refer to the last digit (s). 



Inputs 





6.1 


5.9 


5.7 


Volume, iVf x N T 


24 3 x 48 


16 3 x 32 


12 3 x 24 


Configurations 


200 


350 


300 




1.46 


1.50 


1.57 


K b , m (GeV) 


0.0990, 4.31 


0.0930, 3.73 


0.0890, 2.87 


K c , m Q (GeV) 


0.1260, 1.07 


0.1227, 1.05 


0.1190, 0.96 


K a , mo (GeV J 


0.1373, 0.092 


0.1385, 0.091 


0.1405, 0.093 


K„ 777 n ( CtP\M 
r^q 5 "HI V / 




1382 107 
0.1388, 0.075 
1391 059 
0.1394, 0.043 


0.1410, 0.076 
1415 059 
0.1419, 0.045 




Elementary output 


s 




^crit 


0.13847+2 


0. 14021 ± 3 


0.14327^2 


<hp-is (GeV) 


2.641S 


1.801 6 7 


1.16+J 


a- f l (GeV) 


2.40+™ 


1.471 J 


0.89+2 


2tt/N S (l (GeV) 


0.686 


0.707 


0.607 


Uq 


0.8816 


0.8734 


0.8608 


ocv{2/a) 


0.171 


0.192 


0.227 



III. LATTICE CALCULATION 



This work uses three ensembles of lattice gauge field configurations, which have been 

and 



used in previous work on heavy-light decay constants [ 28|, 14 1 , light-quark masses [2S 



quarkonia BC 



The quark propagators are the same as in Ref. but we now use 200 
instead of 100 configurations on the finest lattice (with f3 = 6.1). The input parameters for 
these fields are in Table |, together with some elementary output parameters. 

The quark propagators are computed from the Sheikholeslami-Wohlert action, which 
includes a dimension-five interaction with coupling csw- F° r heavy and light quarks we 
adjust csw t° the value Uq 3 suggested by tadpole- improved, tree-level perturbation theory, 
and the so-called mean link uq is calculated from the plaquette. The hopping parameter k 
is related to the bare quark mass. For bottom and charmed quarks, Kb and k c are adjusted 
so that the spin-averaged kinetic mass of the corresponding IS quarkonium states match 
experimental measurements. For light quarks, k s and K q are fixed from light meson spec- 
troscopy, using leading-order chiral perturbation theory and the experimental kaon and pion 
masses. We also list the tadpole- improved bare quark mass in GeV, 



(3.1) 



mod = — — 



Uq \2k 2/? cr it 



where the critical quark hopping parameter K crit makes the pion massless. Although this 
mass is just a bare mass, it shows that the heavy quarks are heavy, and the light quarks light. 
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We calculate the three-point function in Eq. fl2.20|) with degenerate spectator and daugh- 
ter light quarks. At each lattice spacing we have propagators corresponding to the strange 
quark. We refer to this decay as B s — > r] s lu, writing r] s for the pseudoscalar ss state in 
analogy with quarkonium. At (3 = 5.9 and 5.7 we have additional light quark propagators, 
with hopping parameter n q , covering the range \m s < m q < m s . 

The lattice spacing a in physical units must be set through some fiducial observable. 
As a rule [30] we prefer the spin- averaged IP- IS splitting of charmonium, Amip_is- For 
comparison we give the value of a' 1 defined through the pion decay constant f n . The 
discrepancy means that Am 1P _i S / f n does not agree with experiment; this is thought to be 
largely due to the quenched approximation, because it remains even as a is decreased. 



The renormalized strong coupling «y(2/a) at scale 2/ a is determined as in Ref. |27]]. It 
is an ingredient in the calculation of the short- distance coefficients py.. and pv ± , introduced 
in Eqs. (|2T7D and ( gl8|) . 

In the three-point functions the heavy-light meson is at rest, while the momentum of the 
light daughter meson is varied. In a finite volume only discrete values of spatial momentum 
are accessible. We compute the three-point function with = 2-n-n/Nsa, for integer mo- 
mentum n G {(0, 0, 0), (1, 0, 0), (1, 1, 0), (1, 1, 1), (2, 0, 0)}. As one can see in Table Q one 
unit of momentum is about 0.7 GeV in the boxes used here, so our calculations cover the 
range < p < 1.5 GeV. 

We obtain the energies, matrix elements, and Zh factors by fitting Eqs. ( |2.20|) and (|2.22|) 
with a ^-minimization algorithm. Statistical errors, including the full correlation matrix 
in x 2 , are determined from 1000 bootstrap samples for each best fit. The bootstrap procedure 
is repeated with the same sequence for all quark mass combinations and momenta, and in 
this way the fully correlated statistical errors are propagated through later stages of the 
analysis. 

The right-hand sides of Eqs. (|2.20|) and (|2.22|) are the first term in a series, with another 
term for each radial excitation. We reduce contamination from these states two ways. First, 
we keep the three points of the three-point function well separated in (Euclidean) time. The 
light meson creation operator O v is always at t i — and the heavy-light meson annihilation 
operator at tf = Nt/2. We then vary the time t s of the current and the range At of time- 
slices kept in the fit, to see when the lowest-lying states dominate. The final choice is made 
by demanding that x 2 /d.o.f. is acceptable and, then, minimizing the statistical errors while 
maximizing At. For acceptable fits we have 3 < At < 6. The extraction of the desired 
matrix elements is shown in Fig. ^] for several light-meson momenta and typical quark mass. 
The best fit and error envelope are indicated by the solid and dotted lines respectively. The 
second way to isolate the lowest-lying states is to choose interpolating operators, Ob and O n 
in Eq. (|2.19|) , to have a large overlap with the desired state. This is done by smearing out 
the quark and anti-quark with IS and 2S Coulomb-gauge wave functions, as in Ref. We 
also examine point-like, or 5 function, operators, but for light mesons at higher momenta 
we find that the 5 source does not yield good plateaus [02][ . The different combinations of 
sources and sinks allow us to check explicitly for excited state contributions by comparing 
results from fits with different smearing functions. Figure [3] compares results for the matrix 
element (B s \V^\r) s ) at n — (1,1,0), obtained from IS source and sink and from IS source 
for the light meson and 5 sink for the B s . The 1S-1S correlation functions yield the cleanest 
matrix elements, so we take our central values from them. 
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FIG. 3. Isolation of the lowest-lying states with different smearing functions, for n = (1, 1,0) 
and quark masses as in Fig. ^. The solid symbols have the standard IS source and IS sink; the 
open symbols have a 5- function sink for the B s . 

IV. ANALYSIS OF FORM FACTORS 



From the exponential fits to three-point correlation functions described in Sec. [TTJ we have 



the matrix element, (7] q (p)\V tl \B q (0)) , for quark masses m q ^ m s and final-state momenta 
\p\ < 1.4 GeV. We must now extend these data to lower quark mass, until the mass of 
the qq pseudoscalar reaches the pion mass. Furthermore, the more important form factor 
f±(E), which is essentially (r) g (p) \ Vj\B g (0)) /pj, is directly calculated only for non-zero three- 
momentum. In the finite volume used here, the lowest non-zero momentum is already 
0.7 GeV, and we would like to extend to lower values, calling for another extrapolation. 

The extrapolation in quark mass can be guided by chiral perturbation theory. To ex- 
trapolate in momentum, however, there is no firm theoretical guide, so we must exercise 
caution. Fortunately, this extrapolation is problematic only in the kinematic regime where 
phase space suppresses the rate. Consequently, neither extrapolation introduces a model. 
We also have checked that the order in which the momentum and chiral extrapolations are 
done has no significant effect the final result. 



A. Momentum interpolation and extrapolation 

Ultimately, we want to compare results at the three different lattice spacings. There- 
fore, we interpolate the lattice data to a fixed set of physical momenta. To start, we 
convert the lattice data to physical units using a!p_ is . Figure f| shows the underlying 
data for B s — > r\ s lv at (3 = 5.9 and 6.1, along with interpolated points. The verti- 
cal (horizontal) error bars on the underlying data come from the statistical uncertainty 
in F\\ ; ± (a -1 ). We interpolate log/j_a -1 / 2 (log /||a +1 / 2 ) linearly (quadratically) in p 2 to 
\p\ e' {0,0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0, 1.1} GeV. This set forms the basis of all 
further analysis. The statistical error bars of the interpolated points are vertical only, be- 
cause both statistical errors are propagated through the interpolation. 
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We must extend the interpolation to an extrapolation to obtain an estimate of f± for 
p < 0.7 GeV. As the pion becomes softer and lighter one expects from Eq. ( |2.5|) that the 
dependence on E (and hence p) is sensitive to the Ansatz for extrapolation. The B* pole 
gives f± a peak at low momentum, and the height of the peak rises as the quark mass 
decreases. This shape is hard to capture, as is shown in Fig. |5|, unless the fit is constrained 
to it. For p > 0.7 GeV the pole fit agrees perfectly with the method described above. But 
as p is decreased into the region of extrapolation, the two forms start to deviate. Above 
0.4 GeV the agreement is still good, so we make a cut here. For smaller momenta phase 
space suppresses the number of events, so this cut has no serious ramifications. For D 
decays the situation is much the same, as shown in Fig. || Therefore, we impose the same 
low-momentum cut in this case. Other functional forms, such as rational, do not make much 
difference in dT /dp cx p A \f + \ 2 /E, once the cut at p m \ n = 0.4 GeV is imposed. 

At high momentum there are other difficulties. The signal-to-noise ratio of the three- 
point function deteriorates. For the highest momentum, n = (2,0,0), we cannot always 
extract a convincing matrix element: in some cases the plateau At is just 2 time-slices, and 
three-point functions with different sources and sinks do not yield the same value for the 
matrix element. We cannot include these data in the interpolation. For the second-highest 
momentum, n = (1, 1, 1), we cannot extract the matrix elements at lighter m q , so statistical 
errors blow up in the chiral extrapolation. We therefore place a cut at n — (1, 1, 0), which 
corresponds to p max = 1.0 GeV. Indeed, our uncertainties would be smaller with a lower 
upper cut, at the cost of reducing the overlap with the experimental data further still. 



B. Chiral extrapolation 

Following the momentum interpolation, the form factors f\\ and f± at j3 = 5.7 and 5.9 are 
extrapolated to the chiral limit at fixed momentum, guided by chiral perturbation theory. 
From Eqs. (|2.4j ) and (|2.5| ) one can see that the chiral behavior of f» and f± should be very 
different. In particular, /n does not contain a B* pole, at least not at the leading order in 
the chiral expansion. In the form factors, dependence on the light pseudoscalar mass enters 
both through m 2 and E. With our momentum cut, p > 0.4 GeV, and our light meson 
masses, 0.45 GeV < m n < 0.74 GeV, the dependence of E on p remains smooth, so we try 
fits of the form 

f l± = A + Bm + Cm 2 , (4.1) 

where m = log(l + mod). We compare quadratic fits with floating C to linear ones with 
fixed C — 0. The difference in the chiral limit of these different fits is the origin of our 
greatest systematic uncertainty. 

It would be desirable to have quark propagators at lighter quark masses to achieve 
better control on the chiral extrapolation. The computer time would increase substantially, 
however, and the obstacle of exceptional configurations would have to be overcome, for 



example as in Ref. [33 



We note that when p = (or p <C m n ) it would be better [33] to carry out the chiral 



extrapolation at fixed E, instead of fixed p. With p > 0.4 GeV, however, the fixed E 
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FIG. 5. Momentum interpolation (and extrapolation) at j3 = 5.9, for B q — > r/gZz/ and (a) the 
heaviest of our light quarks, with n q = 0.1382; (b) the lightest of the light quarks, with K q = 0.1394. 
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TABLE II. Budget of statistical and systematic uncertainties in this work for the quantities 
T B (0A GeV, 1.0 GeV), T D (0A GeV,0.93 GeV), and T B (0A GeV,0.9 GeV)/T D (0.4 GeV,0.9 GeV). 
All entries in percent. 





To 
-L B 


\V i.\ 

I v ub\ 


± D 


\V j| 

I v cd\ 


To /TV> 




statistical 


+27 

- 9 


+14 
- 5 


+17 

- 8 


+ 9 
- 4 


+10 

- 4 


+ 5 

- 2 


excited states 


D 


q 


D 


q 





q 


p extrapolation 


10 


5 


9 


5 


9 


5 


m q extrapolation 


+16 

-22 


+ 8 
-11 


+ 3 
-18 


+ 2 
- 9 


+13 
- 4 


+ 7 

- 2 


adjusting rag 


6 


3 


2 


1 


8 


4 


HQET matching 


10 


5 


10 


5 


10 


5 


a dependence 


+16 

- 3 


+ 8 
- 2 


+23 
- 6 


+n 

- 3 


5 


3 


definition of a 


11 


6 


4 


2 


8 


4 


total systematic 


30 


15 


+28 
-24 


+14 
-12 


+23 
-20 


+12 
-10 


total (stat © syst) 


+40 
-31 


+20 
-16 


+32 
-26 


+16 
-13 


+25 
-20 


+13 
-10 



extrapolation is probably not essential, although it may reduce the uncertainty from the 
chiral extrapolation. We shall investigate this issue elsewhere. 



V. SYSTEMATIC ERRORS 

As discussed in the previous section, we do not have useful results outside the range 

0.4 GeV <p < 1.0 GeV, (5.1) 

where p = \p n \ is the pion's three-momentum in the rest frame of the B or D. Matrix 
elements with higher momentum are not estimated reliably, and at lower momentum the 
chiral extrapolation used is no longer good. In this section we analyze the systematic 
uncertainties quantitatively, focusing on the partially integrated rates T B (0A GeV, 1.0 GeV) 
and T D (0.4 GeV, 0.925 GeV), defined in Eq. (|Q|), and the CKM matrix obtained from 
Eqs. (|1.4j ) and (|1.5| ). A summary of this analysis is given in Table ||. 



The statistical error is estimated with the bootstrap method, drawing 1000 samples 
for each fit. The bootstrap propagates the statistical uncertainty, including correlations, 
through the interpolation in light meson momentum and extrapolation in light-quark mass, 
so in the end statistics remain a quantitatively important source of uncertainty. 



A. Excited states 

As explained in Sec. |TT|, we take care to isolate the desired lowest-lying n and B states 
from their radial excitations when computing the three-point function of Eq. (|2.19| ). The 
associated uncertainty on the matrix elements (and, thus, the form factors) is computed by 
comparing fits with different smeared and unsmeared interpolating operators. After choosing 
the optimal fit range for each combination of smearing functions, we find deviations in F±_ 
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and F\\ of 1-3 percent, where the high end of the range is for momenta near the upper cut. 
We assign an uncertainty of 6 percent to T B and T D . 

Although we calculate similar matrix elements for each p and for B and D decays, the 
range At of time-slices kept in the fit was chosen independently for each case. Therefore, 
the excited state contamination in T b /Td is partly, but not fully, correlated. A conservative 
error estimate is again 6 percent. 



B. Momentum and chiral extrapolations 

The form factor /+ that enters into the partial width is more sensitive to f± than to f\\. 
Thus, it could be sensitive at small p to the extrapolation described in Sec. |IV Aj . The rate, 
however, is much less sensitive, because phase space suppresses it at small p. For T B the 
variation between linear, rational, and pole forms is ±10%. 

The chiral extrapolation is a major source of uncertainty. Figure |7| shows the chiral 
extrapolation at (5 = 5.9 for fa and f± at n = (1,0,0). We compare three different fits: 
(1) a quadratic fit to the four lightest quark masses; (2) a linear fit to the four lightest 
quark masses; and (3) a quadratic fit to all five light quark masses. The first has the lowest 
X 2 /d.o.f., but the other two are perfectly acceptable. For other momenta the behavior is the 
same. Because the extrapolated result from the first (and best) fit lies between the other 
two, we use it to give our central value, and use the other two as estimates of the systematic 
error. The ambiguity of the fits, and hence the systematic error, could be reduced with 
explicit calculation at smaller m q , but a suitable point is not feasible with our computer 
resources. We are left with an uncertainty of tH% in T B and tj/'o in T D . 

The error bars on the extrapolated points in Fig. |7] show how the statistical uncertainties 
are inflated by the chiral extrapolation. This part of the uncertainty is statistical in nature, 
so it is incorporated into the first line of Table 0. Indeed, it is the main reason the statistical 
uncertainty in T B (T D ) grows from 6 percent (7 percent) with m q = m s to 18 percent 
(13 percent) with m q = m d . 



C. Heavy quark mass dependence 

To examine the dependence on the heavy quark mass we use form factors with a light 
strange quark, because then statistical errors do not mask the effect. Figure |] compares the 
form factors B s and D s decays. There is a significant difference. The quarkonium spectrum 



tunes the (bare) heavy quark mass within a precision of 1-2% [14]], which clearly would have 
no significant effect on the form factors. But because of lattice artifacts in the quarkonium 
binding energy |HJ and because of quenching, the heavy-light spectrum yields a different 



adjustment of bare quark masses. The shift is to lower 1/mpg in Fig. ^. From Eq. 
one sees that f± dominates in / + for B decay. Thus, /+ is smaller with the heavy-light 
adjustment of the bottom quark mass, and T B is 6 percent smaller. On the other hand, f± 
and /|| make a comparable contributions to / + for D decay. It turns out that / + is larger 
with the heavy-light adjustment of the charmed quark mass, and Tp is 2 percent larger. 
The ratio T B /T^ is 8 percent smaller. 
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FIG. 8. The heavy quark mass dependence at several momenta, at j3 = 5.9 for the temporal 
matrix element. The light quarks are strange quarks. The dotted lines are to guide the eye. 
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D. Matching 



As explained in Sec. |T[ our treatment of the heavy quark matches lattice gauge theory 
with Wilson fermions to HQET. This requires calculations of the short-distance coefficients: 
\jvTL2 and l/mg in the effective action; y/ZyuuZybb in the definition of the current; and py.., 
pv ± , and I/TO3 in the description of the currents. As discussed in the previous subsection, 
7712 is adjusted non-perturbatively, by tuning the quarkonium spectrum to agree with exper- 
iment. The normalization factors Zyuu and Z v tb are also computed non-perturbatively, by 
requiring that flavor- conserving matrix elements 7r — > 7r and B — > B, computed by analogy 
with Eq. (|2.2Up , give unit charge. The uncertainty from it is purely statistical and much 
smaller than all other statistical uncertainties. 

The significant systematic effects in the matching procedure come from computing pvn 
and pv ± , and from the mismatch between Eqs. ( |2.9| ) and Eqs. (|2.10|) at the level of dimension- 
four and higher currents. In the present work we do this part of the matching with pertur- 
bative QCD, leading to errors of order a 2 ., a s /rriQ, 1/ttiq, respectively. Let us now consider 
these effects in turn. 

Because they are short-distance quantities, the matching factors py, and py x should 
be calculable in perturbation theory. (Note that all effects that make lattice perturbation 
theory less reliable than continuum perturbation theory are absorbed into y/ZyuuZybb.) We 
have calculated them to one loop, so we write 

p v = I + a s (q*)4np v ] (5.2) 

for p v . and py x . We use the Brodsky-Lepage-Mackenzie (BLM) procedure to choose the 
expansion parameter a s (q*) |36|,|27]]. In the scheme in which the Fourier transform of the 
heavy-quark potential reads V(k) = —Cp4:TTa s (k)/k 2 , the BLM scale q* is given through 

*Ji] 

log(g*a) = J\ (5.3) 

where is obtained from py by replacing the gluon propagator D(k) with D(k) logk 2 a 2 . 
The details of these calculations are similar to those described in Ref. |37], and the results 



are listed in Table [II [] ||38|| . The effects are small for B decays and tiny for D decays. This 
can be understood because the pys are ratios of very similar quantities, so there is good 
cancellation. It is therefore plausible that the two-loop contribution is numerically smaller 
by another factor of a s ~ 0.2, and thus completely negligible. 

Next, we must estimate the uncertainty from the mismatch of the I/tuq term in the 
heavy-quark expansion. This contributes an error on either form factor / 

k/mq f ~ a s bi/ mQ (mQa)mg 1 A QCD / (5.4) 
from 1/7713 an d I/Trig contributions, and ^ m gives the deviation of the short- distance 



coefficients for the lattice and continuum theories. (See Refs. [|ll~l , |T2l for further details.) The 
factor b(mQa) is at most of order unity; for our calculations of D-meson matrix elements it 
is of order m c a < 1. Taking a s ~ 0.2 and Aqcd ~ 500 MeV one finds that these errors, in 
either case, are at most a few percent on / or the CKM matrix. 
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TABLE III. Perturbation theory for matching factors py., and py ± . The one-loop terms p v 



and * py are in units of 10 





Pv 


Pv 


6.1 

a s (q*) 


PV 


n [1] 

Pv 


= 

*JXi 
Pv 


5.9 

a s (q*) 


PV 


Pv 


P = 

*„[1] 
Pv 


5.7 

a s {q*) 


PV 


b 

f\\ 
f± 


536 
-1987 


980 
-3312 


0.159 
0.163 


1.011 
0.959 


817 
-2096 


1591 
-3534 


0.173 
0.181 


1.018 
0.952 


1065 
-2146 


2199 
-3621 


0.196 
0.212 


1.026 
0.943 


c 

f\\ 
f± 


-59 
-947 


13 

-1368 


0.233 
0.169 


0.998 
0.980 


-28 
-1223 


40 
-1821 


0.402 
0.188 


0.999 
0.971 


+63 
-1508 


152 
-2339 


0.184 
0.218 


1.001 
0.959 



Finally, we must estimate the uncertainty from the mismatch of the I/tuq terms: 

*L/m»/ ~ W/m%{mQ a ) m ~Q^Q,Cvf- ( 5 - 5 ) 

There are many contribution at order l/win in the heavy-quark expansion, most of which 
come from iteration of the 1/ttiq terms. Only genuine 1/rriQ terms in the effective action 
and currents can be as inaccurate as Eq. ( |5.5[ ) suggests. Since AQ CD /mg « Aq CD a ~ ot a for 
our lattice data, the error 5i/ m » / is similar in magnitude to that of 8x/ mq f . 

The estimates in Table || derived from Eqs. (|5.4|) and (|5.5|) are very conservative. It 
is plausible that the denominator of heavy-quark expansion is 2m,Q, and it is possible that 
the unknown coefficients are fractions instead of 1-2 as used above. Thus, the matching 
uncertainties may already be negligible. 

The masses of the b and c quarks differ by about a factor of three. The short- distance 
coefficients are functions of rriQa []Tl~| , |l2"| , so the matching uncertainties do not cancel com- 
pletely in the ratio T^/Tp. In particular, on our lattices the mismatch coefficients &i/ m « 
are of order 1 for b quarks, but b\/ mQ ~ m c a and &i/ m a ~ (fn c a) 2 , a s m c a for c quarks. 
Nevertheless, the effects often have the same sign, so we take the uncertainty in the ratio to 
be the same as in numerator or denominator. 



E. Lattice spacing dependence 

For the artificial decays B s , D s — > r\ s lv we have results at three lattice spacings, so 
we can examine how severely the form factors are affected. These decays are good for 
studying the a dependence, because their form factors have small statistical errors. After 
chiral extrapolation, on the other hand, the larger statistical error bars would mask lattice 



spacing effects. Previous experience with decay constants |L4] leads us to believe this will not 
change very much after chiral extrapolation. With the action used in this work the lattice 
spacing dependence is a combination of 0(a s a) and 0(a 2 ) effects from the light quarks and 
gluons, and the a dependence of the heavy-quark short-distance coefficients, discussed in 
the previous subsection. In particular, when the r] s has non-zero recoil momentum p, the 
light-quark lattice effects are 0(a s pa) and 0(p 2 a 2 ). 
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The a dependence of the form factors is shown in Fig. [|. The variation with a is several 
percent, which is comparable to the statistical uncertainty and also to the errors from the 
mismatch of the heavy quark. The observed a dependence is therefore a combination of 
(uncorrelated) statistical fluctuations, lattice artifacts from the light degrees of freedom, 
and from the lattice artifacts [described in Eqs. ( |5.4| ) and ( |5.5|) ] of the heavy quark. They 
cannot be disentangled with the current set of calculations, so it does not make sense to 
extrapolate a — > 0. 

Instead we choose the results from j3 = 5.9, where we have the widest range of light quark 
masses, for our central value and use the other two lattices to estimate the uncertainty. 
Figure [XOj shows the combination p 4 \f + \ 2 /E, which is proportional to dT/dp, at all three 
lattice spacings for B s — > r\ s lv and D s — > T) s lv. As in Fig. | one sees that the dependence 
on a is several percent, and increases with increasing p. By integrating over p we find a 
variation of of _ |% in T B and l 2 g% in T D . 



F. Definition of a and quenching 

Changes in the final results from changing the definition of a can be thought of as a 
crude way to estimate effects of the quenched approximation. In lattice units we obtain 
f\\a l l 2 and /±a~ 1//2 , so converting to physical units introduces a mild explicit dependence on 
the value chosen for a. There is also an implicit dependence that enters through functional 
dependence on E (or p). These two effects are illustrated in Fig. |]. The solid (open) points 
are obtained by defining a so that the IP-IS splitting of charmonium (pion decay constant) 
takes its physical value. The central values in the paper are computed with the IP-IS 
definition. At j3 — 5.9 we repeat the full analysis with the f n definition. We find that T B 
(To) increases by 11 percent (4 percent), and the ratio Tb/Td increases by 8 percent. 

A more serious estimate of the effect of the quenched approximation is impossible without 
generating gauge fields with dynamical quark loops. This would require more computer 
resources than we have at our disposal, and no other group has yet studied these semileptonic 
decays with dynamical quarks. There are results with two light, dynamical flavors for the 
leptonic decay constants /# and fjj, using either lattice NRQCD |39| or our method |T6|J40 



for the heavy quark. In that case one finds an increase of between 10-11 percent over the 
quenched result. 

The exercise of changing the definition of a easily could underestimate the effects of 
quenching. At the same time, we do not expect form factors to be more sensitive than 
Thus, a provisional estimate of a uncertainty in Tb,d of 10-20% seems reasonable. 



VI. RESULTS 

The main results of this paper, given in Eqs. ( |1.6| ) and ( |1.7| ), are the quantities 
T B (0A GeV, 1.0 GeV) and T D (0A GeV, 0.925 GeV), which are proportional to the par- 
tially integrated rates. It may also be of interest to present the results in other ways. In 
this section we give results for the ratio Tb/Td, as well as results for Tb, Td, and Tb/Td 
with a lower upper cut. We also give results for the form factors themselves. 
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FIG. 10. Comparison of p A \f + \ 2 /E at three lattice spacings for (a) B 
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TABLE IV. Budget of statistical and systematic uncertainties in this work for the quantities 



Tb (0.4 GeV,0.8 GeV), 


To (0.4 GeV,0.< 


3 GeV), 


and T B (0.4 GeV, 0.8 


GeV)/T D (0.4 GeV,0. 


8 GeV). 


All entries in percent. 














uncertainty 


Tr 


V„h\ 

v UO \ 


Tn 
1 D 


| 'col 


Tr/Td 1 


Vn.h 1 V nA 
v UO 1 ' CCL \ 


statistical 


+21 
-10 


+n 

- 5 


+18 

— y 


+ 9 

r 

— O 


+10 

- 5 


+ 5 

- 3 


excited states 

V..'-/\_ V. ± \ V VI > ' U C. i, LVU 


4 


2 


4 


2 


4 


2 


p extrapolation 


8 


4 


10 


5 


6 


3 


nriq extrapolation 


+14 
—22 


+ 7 

— 11 


+11 

-11 


+ 6 
_ g 


+13 

— 8 


+ 7 

— 4 


adjusting mn 


3 


1 


4 


2 


6 


3 


HQET matching 


10 


5 


10 


5 


10 


5 


o dependence 


+14 
- 2 


+ 7 
- 1 


+17 

- 8 


+ 9 
- 4 


5 


3 


definition of a 


9 


5 


3 


2 


8 


4 


total systematic 


+26 
-28 


+13 
-14 


+26 
-21 


+13 
-10 


+21 
-19 


+n 

- 9 


total (stat © syst) 


+33 
-29 


+17 
-15 


+31 
-23 


+16 
-11 


+23 
-20 


+12 

- 9 



Many uncertainties cancel in the ratio of B and D rates: the statistical error is correlated, 
and the systematic errors are similar in nature. Because of heavy-quark symmetry it is most 
sensible to form a ratio with the same cuts for both. We find 

T B (0.4GeV,0.9GeV) _ + . 20 +0 . 26 



T D (0.4 GeV, 0.9 GeV) 



2.04t ° ™ tom ± 0.10 ± 0.20 ± 0.29, (6.1) 



where the uncertainties are from statistics, chiral extrapolation, a dependence, HQET 
matching, and other miscellaneous sources. A more detailed budget of the last uncertainty 
is given m Table 0. 

As mentioned above, raising the upper cut p max increases the uncertainty. Conversely, 
lowering p max decreases the uncertainty. Repeating the full analysis at p max = 0.8 GeV, we 
find 

T B {0A GeV, 0.8 GeV) = 0.294+°;^ ±g;gg ± 0.029 ± 0.038 GeV 4 , (6.2) 
T D (0A GeV, 0.8 GeV) = 0.1451 ^ ± 0.0161 ^ ± 0.014 ± 0.017 GeV 4 . (6.3) 



and the ratio 

T B (0.4 GeV, 0.8 GeV) n ,,- 

- 2.03lg;Jg t°of e ± 0.10 ±0.20 ±0.24. (6.4) 



T D (0.4 GeV, 0.8 GeV) 



Table |IV| shows a budget of systematic errors, similar to Table ||. As one can see from 



comparing the last two lines in Tables [II] and [IV|, the total uncertainty is several percent 
lower with p max = 0.8 GeV. 

Finally, we give our results for the form factors. Table gives the results for form 
factors in the decay B — > nlu. Listed are f\\ and fx, which emerge directly from our lattice 
calculations, and / + and /o, which appear in the expression for the differential rate. In 
every case the first error is statistical and the second adds the systematic uncertainties in 
quadrature. Table lists the same information for D — > ixlv. Our final results for /y and 
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TABLE V. Form factors with statistical and total systematic errors for the decay B — > irlv. 
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TABLE VI. Form factors with statistical and total systematic errors for the decay D — > wlv. 
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f± are obtained at (3 = 5.9, after chiral extrapolation, with the systematic errors estimated 
as described in Sec. 0. In particular, the estimate of lattice spacing effects uses results from 
all three lattice spacings. For p < 0.4 GeV our extrapolation of f± in the pion momentum 
is no longer reliable, so we do not quote results for it. 

The physical form factors f+(E) and fo(E) are obtained from Eqs. fl2.2|) and (|2.3|) 
using the tabulated results for f±(E) and f\\(E), physical meson masses, and energy 



E = ym 2 +p 2 } for each p = \p w \ in our set of pion three-momenta. They are shown 
in Fig. |TT[ Tables [V] and [VT] also include the combination p A \f + \ 2 /E; for massless final-state 



leptons the differential rates are given by 

= \V ub \ 2 (2.9328 ps^GeV" 4 ) (6.5) 
dp E 

= \V cd \ 2 (1.0358 ps^GeV" 4 ) t\l±L. ( 6 . 6 ) 
dp E 

Other differential distributions can be obtained from the latter by changing variables with 
dp/dE = E/p and dp/dq 2 = E/2rriBP or E/2mr>p. From Eqs. (|6.5|) and (|6.6|) one sees that 
the phase-space factor p 4 suppresses the rate in the low-momentum region where we cannot 
quote /+. 



VII. COMPARISON WITH OTHER RESULTS 



In this section we compare our results to recent published |19,20| and preliminary [21 



work from lattice QCD. The comparison is apt, because three different methods for treat- 
ing the heavy quark on the lattice have been employed. We use Wilson fermions with the 



SW action, normalized to have a consistent heavy-quark limit. References ||19|j20|] use Wil- 
son fermions (with the SW action and light-quark normalization conditions) at nig near 
and below the charmed quark mass, and extrapolate up to m&. Reference [21 uses lattice 
NRQCD |Il[] (with the power-counting of HQET [(42]) and, as we do, calculates the form 
factors directly at the bottom quark mass. 



Figure [12] shows results from Refs. fL9| , pO"f together with ours. (We do not include results 
from the JLQCD collaboration because they are still preliminary. We anticipate that 
their systematic uncertainties will be similar to ours. At this stage their statistical uncer- 
tainties seem surprisingly large.) Within the quoted uncertainties there is broad agreement 
among the three calculations. There are, however, three noteworthy differences in the anal- 
ysis of the form factors. These are the lattice spacings at which the calculations are done, 
the procedure for chiral extrapolation, and the treatment of the heavy quark. 

Our results are based on lattice gauge fields at three lattice spacings, given in Table |. 
We find that the lattice-spacing dependence of the form factors is mild (with our treatment 
of the heavy quark), even on a relatively coarse lattice at /3 = 5.7. The results of Refs. JB5,Sl| 
are both based on only one set of lattice gauge fields (at (3 = 6.2), whose spacing is slightly 
finer than any of ours. We believe, therefore, that the lattice spacing effects of the gluons 
and light quarks are not a serious source of error, at the present overall level of accuracy, in 
any of the three works. 
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FIG. 11. Momentum dependence of the form factors with all systematic uncertainties included, 
(a) /+ and (b) /o- Squares (circles) denote B (D) decays. Solid symbols are independent of the 
momentum interpolation. 
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FIG. 12. Comparison of the differential decay rate (at fixed lattice spacing) vs. E: open dia- 
monds [p|, open triangles [ 20 1 , and solid squares (this work). 
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In our work, the largest uncertainty comes from the chiral extrapolation at fixed pion 
momentum p. As explained in Sec. VB , this uncertainty arises because the linear and 
quadratic fits give moderately different results. References |T9| , pO| do not have enough values 
of the light quark mass to be able to check whether a term quadratic in m q is needed to 
describe their data. Because those works extrapolate at fixed q 2 , however, it is plausible 
that the curvature seen in Fig. [7| would go away, and that a linear chiral extrapolation would 
then be adequate. 

The interpolation in pion momentum, or energy, is another difference. It leads to the 
apparent difference in the shape of the spectrum in Fig. |T2|. If we choose the pole form 
suggested (for small E) by Eqs. Q2.4| ) and (|2.5|), the shape of the spectrum is less humped, 
though not as flat as the spectra from Refs. |l9|,20i . In those works a pole Ansatz different 
from Eqs. ( |2T4] ) and ( |2.5| ) was used. 

The most significant difference in the three calculations is the treatment of the heavy 
quark. Although the same action and similar currents are used, the bare quark mass and 
the normalization of the current are adjusted differently. The normalization conditions 
chosen in Refs. |T9| , pO| are designed for the rriQa — > limit, and at finite rriQa they leave 
systematic uncertainties of order (rriQa) 2 . To reduce these, Refs. [fL9| , p0f calculate with the 
heavy- light meson mass near and below 2 GeV and extrapolate up to mg. This procedure 
leads to their largest quoted uncertainty. The statistical error increases, as it must in any 
extrapolation. There are also systematic effects, which are estimated by trying linear and 
quadratic fits in 1/mg. For at least two reasons, this test may underestimate the systematic 
uncertainty of the extrapolation. First, the compatibility of the fits shows only that the 
dependence on 1/mg is smooth in the employed range of the quark masses. It does not 
show that the heavy-quark expansion is reliable below 2 GeV. This problem is especially 
severe for Ref. [T^j, which has heavy-light meson masses as low as 1.2 GeV. Second, the 
lattice artifacts of order (m^a) 2 may well be amplified by the extrapolation. This problem 
would be especially severe for Ref. PU| , which has mga as high as 0.7. 

Our normalization conditions coincide with those above as rriQa — > 0. At finite rriQa, 
however, they are chosen to eliminate lattice artifacts that grow with rriQa [11]. This is made 

as 



possible by using HQET to match the lattice action and current to continuum QCD flT2 
reviewed in Sec. f|. The advantage is that, as with lattice NRQCD f41,21[, the calculations 
can be done directly at uiq = nib, without an extrapolation in 1/rriQ. Of course, we must 
assume that HQET is valid for the b quark, but that is safer than assuming that it is valid 
for m Q ~ 1-2 GeV. 

A feature of our approach is that it leads to a somewhat complicated pattern of heavy- 
quark discretization errors. There is, however, a corresponding pattern of systematic un- 
certainties in the results of Refs. |i9|,|2~0"|j . In particular, there are corrections to the normal- 



ization of order a s (rriQa) and of order {rriQa) Aqcd/^q = rnQdAQCDa in the 1/mg term 
of the heavy-quark expansion. Estimates of the magnitude of these errors — before or after 
1 / mq extrapolation — are absent from Refs. []l9| , p0| . (The corresponding errors in our work, 
which we address quantitatively in Sec. |V PL are of order a 2 s in the normalization and of 
order c^AgcD/"^ in the 1/niQ term.) 

The calculation of semileptonic form factors for B — > tiIu, and related D decays, has 
also been carried out with quark models and QCD sum rules. At present the uncertainties 
from lattice QCD are comparable to those based on light-cone sum rules [|3|,4~3]. The latter 
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TABLE VII. Strategies for reducing statistical and systematic uncertainties of semileptonic 
form factors. 



uncertainty 



strategy 



statistical 
excited states 
p extrapolation 
m q extrapolation 
adjusting rag 
HQET matching 
a dependence 
definition of a 
quenching 



more lattice gauge fields; better m q extrapolation 
longer time extent; better operators 
larger finite volume; better statistics 

better statistics; more values of m q ; fixed-E 1 extrapolation 
unquench 

match up to a s A QC D/"iQ, (Aqcd/"Iq) 2 

more lattices 

unquench 

unquench 



have the advantage that they are most applicable for energetic pions, so (for B decay) they 
overlap better with th distribution of events in an experiment. On the other hand, it seems 
difficult to reduce the uncertainties from sum rules down to the level of a few percent, which 
will be needed to match the precision of the B factories. As discussed in the following 
section, however, all uncertainties of the form factors are reducible with lattice QCD. 



VIII. CONCLUSIONS 



In this paper we have presented results for the form factors and differential decay rates for 
the semileptonic decays B — > ttIu and D — > nlu. The total uncertainties are 30-35 percent 
(for the rate) and, hence, would yield a theoretical uncertainty to the CKM matrix of 15- 
18 percent. We have attempted a complete analysis of the systematic uncertainties, at 
least within the quenched approximation. A rough estimate of the additional error from 
quenching is 10-20 percent (on the rate). 

A more important, thought less specific, result of this paper is a demonstration that, 
within the quenched approximation, all uncertainties are controllable. Table [VI]] gives a 
sketch of what is needed to reduce all sources of uncertainty. In almost every case, the 
remedy is simply more computer time. That, in fact, is promising, since the computer 



used in this work is already ten years old. Given the experience of the CP-PACS |T6| , |39 
and MILC [E0|| collaborations with heavy-light decay constants, it should be feasible to 
repeat our analysis on a modern supercomputer with unquenched gauge fields. In summary, 
there do not appear to be any technical roadblocks to reducing the uncertainties in lattice 
calculations to a few percent or better, over the course of the present round of experiments. 

In the case of the uncertainty labeled "HQET matching" better calculations of the various 
short- distance coefficients introduced in Sec. || will be needed to be sure that the total 
uncertainty is only a few percent. This is not a computational problem but a theoretical 
one, which arises also in calculations with lattice NRQCD. The alternative would be to 
reduce the lattice spacing dramatically, so that rriQa and Aqcd/^q can be simultaneously 
small. But, since computer requirements scale as a -1 to a high power, that would seem 
to be a long way off. One would also have to sacrifice some other improvements, such as 
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removing the quenched approximation. 

For semileptonic decays of B orD to vector mesons such as p and 00 more study is needed. 
The vector mesons decay hadronically, and in the quenched approximation these decays are 
absent. Even in unquenched calculations, however, there are still issues that may need 
explicit analysis. In particular, if the calculations are done at largish light quark masses, 
the decay may be kinematically forbidden. Because there is not yet much experience with 
unquenched calculations, it is not yet clear whether one can smoothly extrapolate vector 
meson properties from this region to the physical meson masses. It is, thus, hard to anticipate 
how well lattice QCD will do here. This is unfortunate, because the experimental errors for 
semileptonic decays into vector mesons are expected to be somewhat smaller. 

In any case, semileptonic decays of B mesons ultimately will provide one of the most 
accurate constraints on the unitarity triangle, through a determination of \V u b\. Indeed, 
if new physics lurks behind B°-B° mixing, it is essential to have constraints on the CKM 
matrix through charged- current interactions like b — > c and b — > u. 
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